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^j ! Abstract 

We present a comprehensive study of the masses of pseudoscalar and vector 
mesons, as well as octet and decuplet baryons computed in 0(a) improved 
quenched lattice QCD. Results have been obtained using the non-perturbative 
definition of the improvement coefficient c sw , and also its estimate in tadpole 
improved perturbation theory. We investigate effects of improvement on the 
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incidence of exceptional configurations, mass splittings and the parameter J. 
By combining the results obtained using non-perturbative and tadpole im- 
provement in a simultaneous continuum extrapolation we can compare our 
spectral data to experiment. We confirm earlier findings by the CP-PACS 
Collaboration that the quenched light hadron spectrum agrees with experi- 
ment at the 10 % level. 
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I. INTRODUCTION 

Despite recent efforts in simulating lattice QCD with dynamical quarks |I|-§[ the 
quenched approximation is still widely used. While precision tests of QCD through nu- 
merical simulations with dynamical quarks are not possible with the present generation of 
machines, accurate calculations of experimentally known quantities, such as the light hadron 
spectrum, can be performed using the quenched approximation. Recently the results of such 
a benchmark calculation using the Wilson fermion action have been presented by the CP- 
PACS Collaboration ji|, superseding a similar study performed earlier by GF11 ||. Results 
from a similar calculation employing staggered fermions were published in |J. In ref. 0] it 
was concluded that the quenched light hadron spectrum deviates significantly from experi- 
ment by about 10%. 

In order to reach this level of precision one needs to have control over many systematic 
effects, in particular lattice artefacts. In refs. 0J^] extrapolations to the continuum limit 
were performed, thus eliminating the dependence on the lattice spacing a. However, since 
the leading cutoff effects for Wilson fermions are linear in a, it is desirable to corroborate 
these findings and extend the analysis to weak hadronic matrix elements by performing a 
similar study using an improved action. 

To leading order in a the Symanzik improvement programme amounts to adding the 
well-known Sheikholeslami-Wohlert term to the fermionic Wilson action |7| 

SS = -c sw — Yl i'(x)a^F^(x)ij(x). (1) 

Provided that c sw is chosen appropriately, spectral quantities such as hadron masses approach 
the continuum limit with a rate proportional to a 2 . Non-perturbative determinations of 
c sw have been performed in the quenched approximation |8],|9|] and for n/ = 2 flavours of 
dynamical quarks [fLOX Estimates of c sw in tadpole improved perturbation theory |Ll] are 
also widely used. Results for quantities in the light hadron sector using one or the other of 
the two methods have appeared recently [fl~2| p~5|,|9|,p~6|^T~^1 . 

In this paper we present results for the quenched light hadron spectrum in the continuum 
limit, using data computed for both non-perturbative and tadpole improved definitions of 
c sw at several values of the lattice spacing. By combining the two datasets and performing 
a simultaneous continuum extrapolation, we obtain an independent check of the results 
reported in [|,[5| , using a different discretisation of the theory. Here we concentrate on the 
light hadron spectrum. Our results for weak matrix elements such as decay constants will 
be published elsewhere. 

The main conclusion of this work is that the previously observed agreement of the 
quenched light hadron spectrum with experiment at the level of 10 %, is confirmed. Further- 
more, we present qualitative and quantitative analyses of the effects of 0(a) improvement on 
mass splittings, the parameter J, the quark mass dependence of hadrons and the approach 
to the continuum limit. In many ways this work is a continuation of a previous paper fl2| . 

The outline of the paper is as follows: in sect. f| we present the details of our simu- 
lations, including the definition of improvement coefficients and our numerical procedures. 
Sections [T (I] and [V] contain discussions of the "raw" results in the mesonic and baryonic 



sectors, respectively. The quark mass dependence of hadron masses is discussed in sect. [V]. 
In sect. |Vl]we present our results extrapolated to the continuum limit. Detailed comparisons 
of our results and conclusions are presented in sect. |V1J| . 

II. DETAILS OF THE SIMULATION 

A. Improvement coefficients and simulation parameters 

We have generated gauge field configurations using the Wilson plaquette action at three 
values of (3 = 6/g 2 , namely j3 = 5.7,6.0 and 6.2. We used the same hybrid over-relaxed 
algorithm described in |Tj|. For the fermions we have used the 0(a) improved Wilson 
action defined by 

s^[u3M = sf[u3M 

- c swy J2 ^(^a^F^ix^ix), (2) 

where S"p is the standard Wilson action and F^ u is a lattice definition of the field strength 
tensor. The improvement coefficient c sw has been calculated to one loop in perturbation 



theory p0,21 



c sw = 1 + 0.267^ + O((? 4 ). (3) 

It has also been determined non-perturbatively for (3 > 6.0 in ref. || and for (3 > 5.7 in j9]. 
We have computed quark propagators at (3 — 6.0 and 6.2, using the non-perturbative 
determination of c sw from || 

l-0.656gg-0.512gg-0.054gg 

1-0.922^ ' ^-°- u - W 

Furthermore we have used tadpole improved tree- level estimates for c sw , 

e d =«o" 3 > ^ = i(ReTrf/ P ), (5) 

in order to calculate quark propagators at (3 — 5.7, 6.0 and 6.2. In the following we shall 
refer to the datasets computed using either c^ or c*^ d as NP and TAD, respectively. 

Our values for the hopping parameter k were chosen such that they straddle the region of 
the strange quark mass. The simulation parameters for each dataset are compiled in Table |, 
which also contains the estimates of the spatial extensions for each lattice in physical units. 
Exceptional configurations which were encountered at /3 — 6.0 for c^ have been removed 
from the statistical ensemble. The incidence of those configurations is examined in more 
detail in subsect. Ill O. 



In ref. |22| it was argued that the bare parameters have to be rescaled in the 0(a) 
improved theory, so that spectral quantities approach the continuum limit with a rate pro- 
portional to a 2 . [] In the quenched approximation the rescaling needs to be performed only 
for the bare (subtracted) quark mass 

mq = ^6"^)' (6) 

where k c is the critical value of the hopping parameter. The rescaled quark mass rh q is 
defined by 

m q = m q (1 + b m am q ) , (7) 

and the improvement coefficient b m has been computed in one-loop perturbation theory as 



i-0MQ2g 2 + O(g*). (8) 



So far b m has been determined non-perturbatively only at f3 — 6.2 |24]]. We have thus used 
the perturbative estimate in eq. (||), evaluated with a "boosted" coupling g 2 = g 2 /uQ unless 
stated otherwise. In practice we found that the details in the evaluation of b m (e.g. bare 
versus boosted perturbation theory) have little influence on our results. 

B. Hadron correlators and fitting procedure 

Our quark propagators were calculated using both local and smeared sources and sinks. 
The smearing was performed using either the "fuzzing" technique described in ref. [2JJ or the 



Jacobi smearing algorithm of ref. [2q|. Both smearing procedures are gauge invariant. They 
also have a number of parameters, which can be tuned in order to optimise the projection 
on a given hadronic state. For Jacobi smearing the projection properties are controlled by 
the parameter k$, which appears in the kernel of the smearing operator, and the number of 
iterations, iVj ac PB]| . Based on our experience we always chose us = 0.25 and used iVj ac to 
control the smearing radius. 

The fuzzing algorithm for hadronic correlators has three tunable parameters, denoted 
by c, Nf z and r. The parameter c is the so-called "link-staple mixing ratio" , which appears 



in the construction of fuzzed spatial links (at fuzzing level n) according to [27| 



Uf (x)=V{ djf- x) (x) + J2 u i n ~ 1] ( x ) u j n ~ 1] (a + k) U { r l)] (x + j)}, (9) 



x The rescaling is required if a mass independent renormalisation scheme is adopted in which all 
renormalisation conditions are imposed at zero quark mass. In order for such a scheme to be 
compatible with 0(a) improvement the renormalisation of the bare parameters cannot be avoided. 



where V denotes the projection back into the group manifold of SU(3). The maximum 
number of fuzzing levels is given by iVf z . Throughout this work we have used c = 2 and 
iVf z = 5. The size of the fuzzed source (sink) is then determined by r, which is simply the 
length of the straight path of fuzzed links emanating from the origin into all (positive and 
negative) spatial directions. 

An extensive investigation into the optimal smearing parameters, using the projection on 
both mesonic and hadronic states, was performed at (5 — 6.0 on 16 3 • 48 for c^ |2Sfl . It was 
found that N- ]ac = 16 turned out to be a compromise between good projection properties and 
acceptable noise levels in all types of correlators. Similarly, the optimal radius for fuzzed 
sources was determined to be r = 6. For different (3- values the radius r was scaled with the 
lattice spacing. The type of smearing and the corresponding values of iVj ac or r are listed in 
Table | for all datasets. 

Quark propagators computed using smeared or local sources/sinks were combined into 
hadron correlators. We always use the generic notation 'S' to denote correlators which 
have been smeared, regardless of whether fuzzing or Jacobi smearing was used to smear the 
sources and/or sinks. By 'L' we denote unsmeared ("local") sources and sinks. For instance, 
meson and baryon correlators which have been smeared at the source but not at the sink 
are both labelled 'SL' in this notation. The generalisation to other combinations of source 
and sink smearing is obvious. 

We have computed meson and baryon correlators for degenerate and non-degenerate 
combinations of quark masses. Meson correlators in the pseudoscalar and vector channels 
were analysed, as well as spin- 1/2 (octet) and spin-3/2 (decuplet) baryons. 

Our meson masses were extracted by performing correlated, simultaneous fits to the 
(LL, SS) or (LL, SL) combination of correlators. In most cases we used a double-cosh formula 
to fit the ground state and the first excitation, requiring the masses in the fit formulae to 
coincide for both the LL and SS (or SL) correlators. At j3 = 6.0 on 32 3 ■ 64 the double-cosh 
fits turned out to be unstable, so that we resorted to single-cosh fits to either the SS or SL 
correlator. For baryons we followed the same strategy, using double-exponential fits, and, 
at (3 — 6.0, 32 3 ■ 64, a single exponential. 

All fitting intervals have been determined by performing a "sliding window" analysis, in 
which we first selected the maximum timeslice, t max , of the fitting interval (usually t max < 
T/2) and then pushed i m i n to its lowest value which was compatible with the requirements 
of low x 2 /dof and overall stability of the fitted masses. 

All statistical errors have been estimated using the bootstrap method with 1000 bootstrap 
samples. More details about our implementation of the method can be found in p9| . 



C. Exceptional configurations 

It has been noted that calculations of fermionic quantities occasionally suffer from ab- 
normally large fluctuations, in particular for small quark masses [J30"| |3ll. These fluctuations 
have been linked to exceptionally small eigenvalues of the Dirac operator, and the gauge 
configurations on which they occur are usually called "exceptional configurations". The 
fraction of such configurations in the total statistical ensemble increases for smaller quark 
mass m q and/or larger values of g\, c sw and the lattice volume ||. 



In our simulations we have encountered exceptional configurations at (3 = 6.0, but not at 
the other two j3- values. In order to compare their incidence for c*^ d and c"p we have analysed 
distributions of observables for /3 = 6.0 on 16 3 • 48, using the smallest quark mass in the 
TAD and NP datasets. The chosen observable was the unsmeared (i.e. LL) pseudoscalar 
correlator at t — T/2. 

To this end we have determined the median x m and the values denoting the upper (x u ) 
and lower (x{) ends of the central 68 %. As a measure for the width one can define the ratio 

•Em 

The distributions are quite similar for the TAD and NP datasets. First, their width is 
comparable, since Ax u /x m ~ 0.65 in both cases. Second, both distributions extend smoothly 
out to about x m + 9Ax M . 

There are, however, differences in the tails of the distributions, i.e. the number of 
values encountered far beyond x m + 9Ax u . In the TAD dataset only one configuration is 
encountered, which produces a value at roughly ?>7Ax u above the median, whereas in the 
NP dataset there are three such configurations with values for the pseudoscalar correlator 
at 44Ax u , 65 Ax u and 360 Ax u above x m . 

We draw two conclusions from this analysis. The fact that the width is comparable 
(i.e. the value of Ax u /x m ) suggests that the typical statistical fluctuations do not increase 
significantly as c sw is increased from its tadpole improved perturbative estimate to the non- 
perturbative value. Second, we have confirmed the increase of the fraction of exceptional 
configurations (i.e. those configurations for which the observable shoots up to values which 
are orders of magnitude above the normal level of fluctuations) for larger c sw . The presence 
of a zero eigenvalue of the Dirac operator at a nearby K-value has also been verified for such 



configurations [32 



We did not make attempts to treat exceptional configurations using, for instance, the 



methods described in |3T|. Instead we have chosen to eliminate them from our statistical 
ensemble. That is, at /3 — 6.0, 16 3 • 48 we have removed the two configurations which 
produced the most extreme values in the distribution of the unsmeared pion propagator 
for CgP, and on which the inversion of the Dirac operator did not converge for some of its 
components. The latter also occurred on another configuration if a fuzzed source was used, 
and that configuration was subsequently removed as well. The total number of exceptional 
configurations which were removed for a particular dataset are shown in brackets in Table |. 
Note that no configurations were eliminated from the TAD datasets. 

For our range of (3- values and the corresponding values of c sw the incidence of exceptional 
configurations is still relatively small. Their fraction in the NP dataset amounts to less than 
1% on 16 3 -48 (3% on 32 3 -64), and after the analysis presented here we do not expect serious 
distortions of the statistical ensembles due to their removal. 

III. RESULTS FOR PSEUDOSCALAR AND VECTOR MASSES 

In Tables ||-0 we present our "raw" results for meson masses in the pseudoscalar and 
vector channels, which were obtained from the fits described in the previous section. The 
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fit ranges were determined independently for the (LL, SL) and (LL, SS) combinations of 
correlators. Both combinations gave consistent results, and in general we quote the result 
from the fit which gave the best value of x 2 /dof. 

In the pseudoscalar channel the fits were very stable under variations of the fitting 
interval. By contrast, the fits in the vector channel could in some cases differ by up to one 
standard deviation if a different fitting interval was selected. We estimate the systematic 
error in the mass of the vector meson arising from choosing alternative fitting ranges to be 
at most as large as the statistical error. This systematic error has not been included in the 
tables. 

A. Finite volume effects 

Based on our results for the NP dataset obtained at /3 = 6.0 on 16 3 • 48 and 32 3 • 64, we 
can make a first estimate of finite size effects in the mesonic sector. In physical units the 
spatial extensions of the two lattices correspond toL~ 1.5 fm and 3.0 fm, respectively. 

In the pseudoscalar channel we find evidence for small but significant finite volume 
effects. On the smaller lattice the values for amps are consistently larger. Furthermore, the 
effect shows a trend to increase as the quark mass gets smaller. Both these observations 
are consistent with the expected qualitative features of finite-size effects. The difference in 
amps determined for the two lattice sizes amounts to 0.6% at the largest and 1.5% at the 
smallest quark mass. At all values of k the deviation is a two-er effect. 

By contrast, no statistically significant finite size effects are observed in the vector chan- 
nel. In fact, the values for amy are slightly higher on the larger lattice. This might be 
attributed to the fact that no estimate for contributions from excited states was available 
for all hadron masses computed on 32 3 • 64, since only single-cosh fits could be performed. 
Indeed, one of the caveats in the analysis of finite volume effects in both channels is the 
fact that the data for the two lattice sizes shown in Tables [TI] and [Til] have not been ob- 
tained using the same fitting procedure. We have, therefore, repeated the analysis for the 
smaller volume, by performing single-cosh fits for appropriately chosen intervals. The results 
are consistent with those shown in the tables. However, in the pseudoscalar channel it is 
also possible to choose small fitting intervals close to t max such that the single-cosh fits on 
L rj 1.5 fm produce smaller values with larger errors, which are both compatible with the 
results in the tables and also with those obtained on the larger volume. We conclude that 
the finite size effects observed in the pseudoscalar channel appear to be genuine, but without 
further investigations one cannot rule out entirely that they have a statistical origin. For 
vector mesons no significant effects are observed. 

B. Vector-pseudoscalar mass splitting 

It has been known for some time that lattice results obtained in the quenched approxima- 
tion fail to reproduce the experimental fact that the vector-pseudoscalar hyperfine splitting 
is constant over a wide range of quark masses, i.e. m v — mj>g ~ 0.55 GeV 2 . Indeed, lattice 
estimates for this quantity are in general much lower when unimproved Wilson fermions are 
employed. 

8 



In order to study the effect of 0(a) improvement on the hyperfine splitting we have 
plotted our results for the NP and TAD datasets in Fig. |l|(a) and (b), respectively. In order 
to display the influence of finite lattice spacing, we have expressed our results in units of the 
hadronic radius r [J33H , using its lattice determination in ref. ||34j| . The figure demonstrates 



that with improvement (either non-perturbative or mean field) the hyperfine splittings show 
much weaker variation over the studied range of quark masses, compared to the unimproved 
case (see, e.g. ref. [H5|). However, the small slope in the data for r^(my — mp s ) as a function 



of (romps) 2 suggests that the experimentally observed very weak dependence on the quark 
mass is not reproduced by the lattice data. 

By comparing the NP and TAD datasets it appears that the dependence on the lattice 
spacing is somewhat smaller for non-perturbative c sw . Of course this needs to be corrobo- 
rated in a real scaling analysis at a fixed value of (romps) 2 . 

At first sight it may seem surprising that the lattice results for the hyperfine splitting 
overestimate the experimentally observed values. However, lattice results for the splittings 
in the quenched approximation in physical units depend strongly on the choice of scale. 
Indeed, if the scale is set using uik the lattice values are much closer to experiment |35],|28| . 
However, the main focus of this discussion is the analysis of the dependence on the quark 
mass and the lattice spacing. 

C. The parameter J 



The parameter J was introduced |36| as a means to detect deviations between the 
quenched approximation and the observed hadron spectrum without relying on chiral ex- 
trapolations. It is defined through 

J = **£?-, m v /m PS = mK ./ mK , (11) 

and is thus related closely to the slope of the vector-pseudoscalar splitting discussed above. 
Its phenomenological value has been determined from the experimentally measured masses 
as J = 0.48(2). 

In Fig. |2] we plot our results for all our datasets as a function of the lattice spacing. Our 
values confirm previous observations that J is underestimated in the quenched approxima- 
tion. In fact, one finds that the low values for J have little to do with lattice artefacts, since 
there is no sign of the data approaching the phenomenological value for J in the continuum 
limit. We conclude that low lattice estimates for J appear to be an intrinsic feature of the 
quenched approximation. 



IV. RESULTS FOR OCTET AND DECUPLET BARYONS 

Our results for masses of octet and decuplet baryons are shown in Tables |VI|- |X11| . They 
have been obtained by performing double-exponential fits to the (LL,SS) or (LL,SL) com- 
bination of baryon correlation functions, except on 32 3 • 64 for the NP dataset at (3 = 6.0, 
where - in analogy to the fits of the large-volume data in the mesonic sector - only single 
exponential fits to smeared correlators were considered. 
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Tables |VT| and |V11| contain the results for baryon masses in the nucleon ( J p = ~ ) and A 
( J p = | ) channels for degenerate combinations of quark masses, together with the fitting 
ranges and the values of x 2 /dof. These results can be extrapolated to the physical values 
of the quark masses in order to determine the masses of the nucleon, the A (1232) and the 
0(1672), as described in sect. |V|. 

A. Baryons for non-degenerate quark masses 

In order to compute the masses of the physical A, E and S states one has to consider 
baryon correlators for non-degenerate combinations of quark masses. In the octet sector 
one has to distinguish between "E-like" and "A-like" correlators. Using the generic notation 
u, d, s to denote quark flavours, we note that E-like states are symmetric in the light flavours 
u, d, whereas A-like states are antisymmetric. On the lattice the corresponding correlators 
are obtained by performing the appropriate contractions |37| , |38| , ^8| . The J p = \ states of 



the E and A are then obtained from the correlation functions by averaging the 11 and 22 
spinor indices. 

The correlators for decuplet baryons, which are symmetric in all flavours u, d, s are 
simpler to construct. They are obtained in terms of the interpolating operator 

D Kijk = e abc (i>?Cj^) n, (12) 

where i, j, k denote the quark flavour, a, b, c are colour indices, and C is the charge conjuga- 
tion matrix. Correlation functions for decuplet baryons are constructed from the correlation 
of D^ij k by projecting out the spin-| component. 

In Tables VTljJ-jXXI we list the results for octet (E, A-like) and decuplet (A-like) baryons 



for non-degenerate combinations of quark masses. The fitting intervals, which are not shown, 
are mostly identical to those chosen for the corresponding channels in the degenerate case 
(Tables |VT] and |V1J| ) . With our statistical accuracy we are not able to distinguish between 
E- and A-like states; the different symmetry properties in the quark flavours corresponding 
to the hopping parameters Ki and k 2 do not manifest themselves in statistically significant 
mass differences. 

By extrapolating or interpolating the data in Ki, k 2 , K3 to the hopping parameters corre- 
sponding to the physical quark masses one obtains the masses of the A, E, E*, S, H* and Q. 
Note that non-degenerate combinations of quark masses have not been computed at j3 = 5.7. 

B. Finite volume effects in the baryonic sector 

The issue of finite-volume effects is of special importance in the baryonic sector where 
these effects are expected to be more severe than for mesons. With our data we can assess 
the influence of finite- volume effects by comparing our results computed at j3 = 6.0 on either 
16 3 • 48 or 32 3 • 64 using the non-perturbative value of c sw . 

The numbers in Tables IV]], [VTlI| and [TX] suggest that on the large volume the mass 
estimates for both octet and decuplet baryons are slightly smaller. For octet baryons this 
decrease amounts to about 1.4% at the largest and 2.7% at the smallest quark masses. The 
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effect is roughly twice as large as for the pseudoscalar mesons discussed earlier. Although 
finite-size effects for octet baryons are not significant at our level of statistical accuracy, this 
does not necessarily indicate that those effects are absent. 

For decuplet baryons the finite-volume effects are more pronounced; they vary be- 
tween 2.4 and 5.5%, again increasing towards smaller quark masses. Here the discrepancy 
between the results on the small and large volumes amounts to about 1.5 standard devia- 
tions. Thus we cannot exclude finite-size effects in our baryon data at a level of up to 2.5% 
for octet and 5.5% for decuplet baryons. 

V. QUARK MASS DEPENDENCE 

In this section we discuss the dependence of mesons and baryons on the quark mass. 
Usually this dependence is modelled using the results of chiral perturbation theory at lowest 
order. It is then quite a delicate problem to decide whether higher orders in the chiral 
expansion have to be included. Furthermore, additional care must be taken in the quenched 
approximation, where one expects deviations from the leading behaviour for very small 
quark masses (i.e. close to the chiral limit), due to the appearance of quenched chiral 
logarithms [^,|y,^,^T[] . We will first motivate the functional forms for the quark mass 
dependence used in this paper, determine the critical hopping parameter and then present 
our results for hadron masses extrapolated or interpolated to the physical quark masses. 

A. Fit ansatz and the critical hopping parameter 

Usually the critical value of the hopping parameter, k c , is determined at the point where 
the mass of the pseudoscalar meson vanishes, mps = 0. The simplest ansatz for the quark 
mass dependence of mps, which is consistent with 0(a) improvement, is 

m PS = B(m qA + m q , 2 ), (13) 

where m qji , i = 1, 2 denotes the rescaled, bare quark mass defined in eq. (|7|). 



Assuming that the ansatz in eq. (0) is justified (i.e. both higher orders in the quark 
mass as well as quenched chiral logarithms are assumed to be absent), we have determined 
n c using both degenerate and non-degenerate combinations of quark masses, by inserting 



the definition of m q into eq. fli"3]), which leads to the general fit ansatz 

™ps = « + P (~ + -) + 7 | 4 + ~ 2 ) , (14) 



where the fit parameters a, f3 and 7 are related to B, k c and b m through 

Bb m 




(15) 

As mentioned in sect. |I| we have used the tadpole improved perturbative estimate at one 
loop for b m . In order to study the sensitivity of k c on b m we have also used its estimate in 
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one-loop perturbation theory in the bare coupling, as well as its tree-level value, 6 m = —1/2, 



and, at j3 = 6.2 for the NP dataset, the non-perturbative determination of ref. ||24j| . In 
order to enable direct comparisons with K c -estimates from earlier simulations using tadpole 
improvement, we have also computed k c for 6 m = for the TAD dataset. Our results, 



which are collected in Table |X11J|, show little dependence on the value of b m . The largest 



deviations are observed at (3 = 5.7. Furthermore, using the non-perturbative value of b m at 
(3 = 6.2 yields a result which is entirely compatible with the ones obtained using one-loop 
perturbative estimates. 

We conclude that for our range of quark masses, estimates for 6 m based on one-loop 
perturbation theory are sufficient to obtain stable results for n c for /3 ^ 6.0. 

We can now justify our ansatz eq. fll3D by plotting mp s as a function of (m q) i + m qj2 )/2 
in units of tq for both NP and TAD datasets. This is shown in Fig. [3] where the lines denote 
the fits based on eq. (|i~3|) . The plots show that no significant departure from the linear 
behaviour predicted by lowest order chiral perturbation theory is observed in the range of 
quark masses investigated. Thus, we find no evidence for higher order terms in the chiral 
expansion of w^ps, nor do our data support the presence of quenched chiral logarithms. The 
latter is most probably due to the fact that the quark masses used in our simulations are 
not light enough. 



Furthermore, we wish to point out that a more sophisticated analysis [|18] of the quark 
mass dependence of the data in Table || revealed that higher-order terms proportional to 
(^q,i — ^q,2) 2 contribute below the 1% level. 

Based on our observations in the pseudoscalar channel, which usually offers the most 
precise information about the quark mass dependence, we have used the following functional 
forms for vector mesons, octet and decuplet baryons: 

m v = Ay + Cy (m q ,i + m q , 2 ) (16) 

m 0ct = A + C (m q ,i + m qi2 + m q , 3 ) (17) 

m Dec = A D + C D (m q ,i + m qi2 + m q , 3 ) . (18) 

The corresponding fits are shown in Figs. ^ and [5] for the NP and TAD datasets, respectively. 
As in the case of pseudoscalar mesons we observe that for our level of precision and range of 
quark masses there is no evidence for curvature in the data. We conclude that eqs. (|T3|-|I"8"|) 
represent appropriate fitting functions for our data. A detailed analysis of more complicated 
models for the quark mass dependence described in p8| resulted in similar findings. 

In a given channel (i.e. pseudoscalar and vector mesons, E-, A- and A-like baryons) 
we have determined the parameters B, Ay, Cy, . . . from uncorrelated, simultaneous fits to 
degenerate and non-degenerate combinations of quark masses. The only exception was the 
dataset at (3 = 5.7, for which only two degenerate combinations of quark masses had been 
computed in the baryonic channels. Therefore, the quark mass dependence at j3 — 5.7 is 
not really controlled. Nevertheless, we have included the results in the following analysis. 

B. Hadrons at physical values of the quark masses 

Our task now is to make contact with the physical hadron spectrum by matching the 
quark masses in eqs. (|T5|-|T%1) to the masses of the physical u, d and s quarks. Here we employ 
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the axial Ward identities, which, for the physical pseudoscalar mesons (in the continuum 
theory), read 

ml± = B{m u + m d ) (19) 

m 2 K ± = B(m u + m s ), m 2 K0 = B(m d + m s ). (20) 

By assuming isospin symmetry, i.e. m„ = m d one can define the so-called "normal" quark 
mass m n through 

m n = \{m u + m d ) (21) 

and the isospin-averaged combination of kaon masses as 

(22) 





ra\ = \{m\ ± + m 2 K0 ), 


such that 






m^± = 2B m n 




m 2 K = B(m n + m s ). 


Inserting m 2 K ± - 


= 493.7 MeV and m K o = 497.7 MeV g2| 



(23) 
(24) 

one obtains^ 

m K = 495.7 MeV. (25) 

Using our lattice estimates for the parameter B we can determine the combination of (bare) 
quark masses m s + m n and the quark mass m n through 

a{m s + m n ) = (aQ) p (26) 

an 

am n = (aQ) n B ^ ' 

where Q is the quantity which sets the lattice scale, and the subscript "phys" denotes that 
the physical ratio is to be taken. It is a well-known fact that in the quenched approximation 
there is an intrinsic ambiguity associated with the lattice scale. In order to estimate this 
ambiguity we have used three different quantities for Q, namely 



Q = r \ 


r = 0.5fm, r^ 1 = 395 MeV, 


Q = m K *, 


m K * =893.9 MeV, 


Q = m N , 


m N = 938 MeV. 



(28) 

Here, the value m K * = 893.9 MeV is the isospin averaged result, i.e. m K * = ^{m K »±+m K *o). 
Lattice data for r /a and its error at all relevant (3- values were taken from ref. | 53fl . In 
particular, we used the interpolating formula, eq. (2.18) of |$4 . 



The following comments apply to our chosen set of lattice scales: 



2 In principle one also has to compensate for the electromagnetic binding energy of about —0.7 MeV 
in eq. (f25j) (see e.g. ref. |1"8||). However, this has not been done in this paper. 
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On the lattice the mass of the K* meson is obtained through an interpolation of 
data points, which is intrinsically a safe procedure. However, the physical K* is an 
unstable hadron with a finite width; resonance effects are not controlled in the lattice 
calculation; 

The nucleon is a stable hadron, but in lattice simulations its mass is obtained only by 
an extrapolation close to the chiral limit. In view of the possible presence of quenched 
chiral logarithms this extrapolation is hard to control. Furthermore, precise lattice 
determinations of baryon masses are more difficult compared to the mesonic sector, 
due to larger statistical errors and the possibility of relatively large finite size effects; 



The hadronic radius tq is known accurately for a wide range of lattice spacings [51 
but a direct experimental measurement is not available. Its phenomenological value 
of r = 0.5 fm is estimated from potential models fitted to experimental data. 



In Table [XIV| we have collected the results for the hopping parameters n n and n s , cor- 
responding to the quark masses am n and am s , obtained for our three different choices of Q 
and using the tadpole improved perturbative estimate for b m in the definition of am^. Thus, 
in spite of the difficulties associated with extrapolations to the chiral limit we have chosen 
to compute am n and to quote lattice estimates for m p , m^ and 1R4. 

The physical vector meson, octet and decuplet baryon masses have been computed by 
inserting the appropriate combinations of am n and am s , corresponding to the physical quark 
content, into eqs. (|T6|-|T8|). The results, in units of either r , mx* or m^ are shown in 
Tables X"V, XVI and |XV1J| , respectively. These tables contain the data with L/a < 24 only. 
For completeness, the non-perturbatively improved data obtained on 32 3 • 64, (3 = 6.0, which 



have not been used in the continuum extrapolation, are listed in Table [XVI I I 



It has been observed that several tests of the quenched approximation can be performed 
at intermediate values of the quark mass (e.g. near the strange quark mass), so that extrap- 
olations to the chiral limit as a reference point are not required [M. Furthermore, it has 
been suggested that properties of the effective chiral Lagrangian can be studied for unphys- 
ical quark masses. A convenient reference point for future lattice studies is then provided 
by the condition (mps^o) 2 — 3.0, which has already been chosen in ref. [|1^]. We have in- 
terpolated our results in the vector meson, E and A channels to that point by defining a 
reference quark mass m re f through 



2(aB) (am rei ) (r Q /a) 2 = 3.0 



(29) 



The results are shown in Table XIX. 



and inserting its value into eqs. 

Another reference point which does not require extrapolations to the chiral limit is defined 
at the point where 

mps/my = 0.7. (30) 

Results for the vector meson and nucleon masses at the reference point defined by eq 



have been quoted in |T4]j9[|. Our estimates in the vector, E and A channels are included 
in Table [X1X| . By comparing the results in Tables |X1X| and [XV| one observes that both 
reference points correspond to the case of degenerate light quarks with masses around that 
of the strange quark. 
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VI. THE CONTINUUM LIMIT 

We are now in a position to discuss the extrapolation of our results to the continuum 
limit. This will finally enable us to make a direct comparison with experimental data and 
the results of ref. ||, obtained using the unimproved Wilson action. 

For hadron masses computed using the non-perturbative determination of c sw the leading 
cutoff effects are expected to be of order a 2 . Detailed scaling studies have confirmed that 
the approach to the continuum limit for spectral quantities P,f43|| and matrix elements |43 



is indeed consistent with such a leading order. By contrast, it is expected that small terms 
of order a cannot be excluded when tadpole improvement is used |TJ . 



From our list of simulation parameters in Table [I] it is clear that separate continuum 
extrapolations of the results for the NP and TAD datasets are not feasible. We have there- 
fore chosen to perform simultaneous extrapolations by assuming leading lattice artefacts of 
order a 2 for the NP and artefacts of both order a and a 2 for the TAD dataset. The ansatz 
for the continuum extrapolation of a generic hadron mass M in units of r$ then reads 

° I n)M| a=0 + A TAD (a/r ) + 5 TAD (a/r ) 2 . (31) 

In other words, one requires that the data computed for non-perturbative and tadpole 
improvement extrapolate to a common continuum value. By combining the results obtained 
at three /3-values for tadpole improvement with those at two /3-values in the case of non- 
perturbative improvement we have five data points to determine the four fit parameters 
M, B NF , A TAD and B™. Note that the data obtained on 32 3 • 64 at /3 = 6.0 have not been 
used. It is worth pointing out that the spatial volume at (3 = 5.7 is larger than those for the 
larger two (3- values which enter the extrapolations. 

In Fig. |^ we show examples of continuum extrapolations based on eq. (|3l|) , namely one 
representative of each of the vector meson, octet and decuplet baryon channels, respectively. 
The extrapolations have been repeated for the other two choices of the lattice scale, i.e. 
Q = rriK* and mjv- The value of x 2 /dof for these fits was quite low (below 1) for all 
channels considered. 

The results are listed in Tables pCV^XVlTI in the row labelled "Cont.". These numbers 
represent the final results for the physical hadrons in this paper. In addition we have also 
performed continuum extrapolations of hadrons interpolated to (mpsro) 2 = 3.0. The results 
are also included in Table [XIX| . 

To check whether or not the continuum extrapolations for different choices of the lattice 
scale are controlled, one can compare, for instance, the continuum result for a particular 
hadron in units of m K *, i.e. (M/ijik*) with the ratio (r M)/(r mK*)- It then turns out 
that not only are the values consistent within errors, but they are also numerically very 
close. This gives us further confidence that the continuum estimates in Tables |XV| - [XVrlj are 
reliable. 

From the results in the tables one can also estimate the size of lattice artefacts at a fixed 
/3-value, say, (3 = 6.0 which roughly corresponds to a ~ 0.1 fm. Using the numbers from 



Table XV one infers that lattice artefacts for both mesons and baryons at /3 — 6.0 are of 



the order of 5 % or less. Apart from Tables |XVHXIX| and Fig. ^ some information about the 
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relative scaling behaviour of the NP and TAD datasets can also be gained from Figs. 0-[5|. 
Here one observes that data for mesons and baryons in units of r are almost independent 
of the lattice spacing for /3 > 6.0. This is particularly pronounced when non-perturbative 
improvement is employed. It is also clear that significant lattice artefacts are present in the 
tadpole improved data at (5 — 5.7. 

VII. DISCUSSION AND CONCLUSIONS 



Our final results in Tables IXVHXVH] and |XIX| can now be compared to other lattice 



calculations and experimental data. 

In ref. |18]| results for vector mesons computed using the non-perturbative value of c s 



were presented. By comparing our results for r my computed at (r mp S ) = 3.0 and those 
for r mx* to ref. |18| we find differences of 1-2 standard deviations at most. Since the 



numerical procedures employed in [113 are quite different to those used in this paper (see 



also |L7|) this agreement is an important check of the stability of our results, both at finite 
lattice spacing and in the continuum limit. 

As mentioned in the introduction, one important goal of this study is to corroborate 
earlier findings by GF11 H and the CP-PACS Collaboration ||, using an improved discreti- 
sation of the QCD action. In Fig. [?| we present our final results in physical units, computed 
using either uik* or m^ to set the scale. Our numbers are compared to the CP-PACS 
results, obtained using m p to set the scale, and to the experimental numbers. 

The first observation is that, on the whole, the two simulations agree quite well, although 
the errors quoted by CP-PACS are in general much smaller. This confirms the previous con- 
clusion that the quenched light hadron spectrum agrees within 10% with experiment. This 
independent confirmation, using results in the 0(a) improved theory, is an important result, 
since lattice artefacts for unimproved Wilson fermions are in general much larger |35||J, so 
that continuum extrapolations can be quite drastic. 

It must be mentioned, though, that the overall precision of our results cannot match 
that of ref. j|] for the following reasons. Firstly, our calculations have been performed in a 
fairly narrow range of quark masses, for relatively small volumes (mostly for L « 1.5 fm). 
Therefore, possible deviations from a linear quark mass dependence close to the chiral limit 
could not be detected or controlled. This affects mainly states like the p, nucleon or the 
A. Indeed, the modelling of the observed downward curvature in the data reported in |l| 
turned out to be a significant factor in the detection of the deviation from the experimentally 
observed spectrum. Furthermore, as far as our baryonic data are concerned, the chiral 
behaviour at (3 = 5.7 is not controlled at all, since only two data points have been computed, 
so that we had to assume linearity in the quark mass. Finally, our data are subject to finite- 
size corrections. For hadron masses computed at the physical values of the quark masses 
these effects amount to at most 2 % for vector mesons, 1 — 3 % for octet and 4 — 8 % for 
decuplet baryons (c.f. Tables |XVHXVHl|) . 



To some extent one can check against possible finite-size effects in our results by setting 
the scale using the nucleon mass. As a baryon the nucleon might also be affected more 
by finite-size effects, and it is reasonable to expect a (partial) cancellation of these effects 
in the ratios mhadron/^iv- Indeed, as can be seen from Fig. |7| our results in physical units 
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decrease when the scale is set by m^ rather than rriK* , although the difference is mostly not 
significant. 

If the differences in the results due to different choices of the lattice scale is not attributed 
to finite-size effects, then they serve to estimate the intrinsic scale ambiguity in the quenched 
approximation in the continuum limit. On the basis of the results in Tables [XV| - [XVII| one 
can infer that in the most extreme cases the difference between the highest and lowest values 
in physical units amounts to about 20 %. We may then assign an uncertainty of ±10 % to 
our results in physical units as a consequence of the scale ambiguity. 

To summarise, we have presented a comprehensive study of the light hadron spectrum 
in quenched QCD, using improved Wilson actions. Qualitative results indicate an improved 
behaviour of the pseudoscalar-vector mass splitting. The parameter J, on the other hand, 
shows no sign of approaching its phenomenological value of 0.48(2), even in the continuum 
limit. This appears to be an intrinsic feature of the quenched approximation. 

Our results show that 0(a) improvement works well for spectral quantities. Hadron 
masses computed for non-perturbative c sw and expressed in units of r show almost no 
dependence on the lattice spacing for (3 > 6.0. We also find that the extrapolations to the 
continuum limit are quite mild in general. That is, for a ~ 0.1 fm lattice artefacts amount 
to about 5 % or less. 

We have presented further evidence that the quenched light hadron spectrum agrees with 
experiment to within 10%. Further technical improvements, including, in particular, the 
modelling of the quark mass dependence and the addition of more points in the continuum 
extrapolations should be implemented to increase the overall precision. 

The next step is the extension of this investigation to decay constants and matrix ele- 



ments. First results have already been presented in |1 1 - J46J . Here, an attractive feature is 
the availability of non-perturbative determinations of some renormalisation factors p7| ) ^5| . 
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TABLES 

TABLE I. Simulation parameters, statistics and smearing parameters for the NP (upper three 
rows) and TAD (lower three rows) datasets. Lattice sizes in physical units are estimated using 
vq to set the scale J34J. The number of exceptional configurations removed from the ensemble is 
denoted in parentheses. 



6.0 
6.2 



L 3 -T 



L [fm] 



#conf. 



smearing 



1.769 16 3 -48 1.5 0.13344,0.13417,0.13455 496(3) 

32 3 • 64 3.0 0.13344, 0.13417, 0.13455 70(2) 

1.614 24 3 -48 1.6 0.13460,0.13510,0.13530 216 



fuzz, r = 6 
jac, iVj ac = 30 
fuzz, r = 8 



5.7 1.568 16 3 -48 2.7 0.13843,0.14077 145 

6.0 1.479 16 3 -48 1.5 0.13700,0.13810,0.13856 499 

6.2 1.442 24 3 -48 1.6 0.13640,0.13710,0.13745 218 



jac, iVj ac = 16 
fuzz, r = 6 
fuzz, r = 8 



TABLE II. Pseudoscalar masses for the non-perturbatively improved data sets. 



p 


L 6 


T 


«l 


«2 


amps [tmin, imax] 


X 2 / dof 


6.0 


16 3 


48 


0.13344 


0.13344 


0.3977±7 3 


6,23] 


23.82 / 30 








0.13417 


0.13344 


0.3553±y 5 


6,23] 


24.18 / 30 








0.13455 


0.13344 


0.3319^ 7 


6,23] 


26.70 / 30 








0.13417 


0.13417 


0.3077^ 8 


6,23] 


26.14 / 30 








0.13455 


0.13417 


0.2805±^ 


6,23] 


27.37 / 30 








0.13455 


0.13455 


0.24931?! 


6,23] 


30.84 / 30 


6.0 


32 3 


64 


0.13344 


0.13344 


0.3952if 


15,31] 


14.49 / 15 








0.13417 


0.13344 


0.3524lio 


15,31] 


14.86 / 15 








0.13455 


0.13344 


0.3284^ 


15,31] 


13.56 / 15 








0.13417 


0.13417 


0.3048±^f 


15,31] 


13.64 / 15 








0.13455 


0.13417 


0.2769ijf 


15,31] 


12.22 / 15 








0.13455 


0.13455 


0.2457t^ 


15,31] 


13.04 / 15 


6.2 


24 3 


48 


0.13460 


0.13460 


0.2803J:io 


8,23] 


30.99 / 26 








0.13510 


0.13460 


0.2492t^ 


8,23] 


29.08 / 26 








0.13530 


0.13460 


0.2361±J| 


8,23] 


28.42 / 26 








0.13510 


0.13510 


0.2149i^ 


8,23] 


31.54 / 26 








0.13530 


0.13510 


0.19981^ 


8,23] 


31.18 / 26 








0.13530 


0.13530 


0.1836^| 


8,23] 


32.04 / 26 
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TABLE III. Vector masses for the non-perturbatively improved data sets. 



p 


L 3 


T 


«i 


K-2 


CLiny l^miri) ImaxJ 


X 2 / dof 


6.0 


16 3 


48 


0.13344 


0.13344 


0.5397l|; 


6,23] 


24.02 / 30 








0.13417 


0.13344 


0.5124l|j 


6,23] 


27.16 / 30 








0.13455 


0.13344 


0.4997+JJ 


6,23] 


29.71 / 30 








0.13417 


0.13417 


0.4852l| 


6,23] 


27.92 / 30 








0.13455 


0.13417 


0.4713l| 


6,23] 


31.96 / 30 








0.13455 


0.13455 


0.4577l| 


6,23] 


30.16 / 30 


6.0 


32 3 


64 


0.13344 


0.13344 


0.5400lg 


10,20] 


13.68 / 9 








0.13417 


0.13344 


0.5143l| 


10,20] 


14.44 / 9 








0.13455 


0.13344 


0.50191^ 


10,20] 


13.75 / 9 








0.13417 


0.13417 


0.48871^ 


10,20] 


13.30 / 9 








0.13455 


0.13417 


0.47621^ 


10,20] 


10.31 / 9 








0.13455 


0.13455 


0.46361^ 


10,20] 


6.69 / 9 


6.2 


24 3 


48 


0.13460 


0.13460 


0.38871^ 


8,23] 


33.55 / 26 








0.13510 


0.13460 


0.37081^ 


8,23] 


28.99 / 26 








0.13530 


0.13460 


0.36451^ 


8,23] 


26.51 / 26 








0.13510 


0.13510 


0.35311^ 


8,23] 


29.56 / 26 








0.13530 


0.13510 


0.347llg 


8,23] 


27.91 / 26 








0.13530 


0.13530 


0.34141^ 


8,23] 


30.98 / 26 



TABLE IV. Pseudoscalar masses for the tadpole improved data sets. 







L 6 T 



«i 



>'<2 



amps 



[train j tmaxj 



x 2 / dot 



5.7 


16 3 


32 


0.13843 


0.13843 


0.73501^ 


6,15] 


23.99 / 14 








0.14077 


0.13843 


0.64041^ 


5,15] 


16.89 / 16 








0.14077 


0.14077 


0.53071^ 


5,15] 


16.00 / 16 


6.0 


16 3 


48 


0.13700 


0.13700 


0.413ll^ 3 


6,23] 


22.57 / 30 








0.13810 


0.13700 


0.3572l^ 6 


6,23] 


21.97 / 30 








0.13856 


0.13700 


0.3320!^ 


6,23] 


27.82 / 30 








0.13810 


0.13810 


0.2927lf 


6,23] 


26.93 / 30 








0.13856 


0.13810 


0.262llf 


6,23] 


29.93 / 30 








0.13856 


0.13856 


0.22681^ 


6,23] 


30.40 / 30 


6.2 


24 3 


48 


0.13640 


0.13640 


0.30331^ 


8,23] 


30.03 / 26 








0.13710 


0.13640 


0.26431^ 


8,23] 


29.29 / 26 








0.13745 


0.13640 


0.2436lJ| 


8,23] 


29.02 / 26 








0.13710 


0.13710 


0.22061^ 


8,23] 


31.97 / 26 








0.13745 


0.13710 


0.1959lie 


8,23] 


30.92 / 26 








0.13745 


0.13745 


0.16801^ 


8,23] 


31.13 / 26 
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TABLE V. Vector masses for the tadpole improved data sets. 



p 


L 3 


T 


«i 


K-2 


CLiny l^miri) ImaxJ 


x7 dof 


5.7 


16 3 


32 


0.13843 


0.13843 


0.9332±$ 


7,15] 


12.80 / 12 








0.14077 


0.13843 


0.86881^ 


7,15] 


11.96 / 12 








0.14077 


0.14077 


0.809±? 3 


7,15] 


10.76 / 12 


6.0 


16 3 


48 


0.13700 


0.13700 


0.53861^ 


7,23] 


23.82 / 28 








0.13810 


0.13700 


0.50301^ 


6,23] 


27.20 / 30 








0.13856 


0.13700 


0.48891^ 


6,23] 


27.83 / 30 








0.13810 


0.13810 


0.4652+^ 


6,23] 


26.92 / 30 








0.13856 


0.13810 


0.4501+^ 


6,23] 


29.23 / 30 








0.13856 


0.13856 


0.4353+^ 


6,23] 


25.43 / 30 


6.2 


24 3 


48 


0.13640 


0.13640 


0.40051^ 


8,23] 


32.65 / 26 








0.13710 


0.13640 


0.37611^ 


8,23] 


28.79 / 26 








0.13745 


0.13640 


0.36481^ 


8,23] 


25.39 / 26 








0.13710 


0.13710 


0.3522i^ 


8,23] 


26.83 / 26 








0.13745 


0.13710 


0.34121^ 


8,23] 


24.85 / 26 








0.13745 


0.13745 


0.33061^ 


8,23] 


28.83 / 26 



TABLE VI. Masses for the nucleoli and A for degenerate quark mass combinations for the 
non-perturbatively improved data. 



p 


L 3 


T 


K 


OWN 


l^min) ^maxj 


x7 dof 


arriA 


[i-min) f-maxj 


X7 dof 


6.0 


16 3 


48 


0.13344 
0.13417 
0.13455 


0.808li, u 

0.711+J! 

0.66511 


[ 9,23] 

[ 9,23] 
[ 9,23] 


25.42/24 
25.74/24 
26.98/24 


0.9131^2 

0.8521:1 
0.7681^ 


[ 9,23] 
[ 9,23] 
[10,23] 


21.84/24 
31.33/24 
40.37/22 


6.0 


32 3 


64 


0.13344 
0.13417 
0.13455 


0.799l^ 
0.700l^ 
0.6411^ 


[ 3,18] 
[ 3,18] 
[ 3,18] 


16.82/12 
15.99/12 
14.63/12 


0.8991^ 
0.818±i| 
0.781±i| 


[ 2,16] 
[ 2,16] 
[ 2,16] 


17.56/11 
20.09/11 
21.66/11 


6.2 


24 3 


48 


0.13460 
0.13510 
0.13530 


0.586l^ 

0.509li[] 

0.487lf 4 


[10,23] 
[10,23] 
[10,23] 


42.46/22 
41.09/22 
30.35/22 


0.6711? 

O.6I8IJ2 
0.596lJ| 


[11,23] 
[11,23] 
[11,23] 


20.35/20 
21.54/20 
20.31/20 
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TABLE VII. Masses for the nucleoli and A for degenerate quark mass combinations for the 
tadpole improved data. 



(3 


L 3 


■T 


K 


am,N 


[£min, '-maxj 


x7 dof 


am/\ 


L^min, tmaxj 


X 2 / dof 


5.7 


16 3 


•32 


0.13843 
0.14077 


1.423lf 
1.183±if 


[ 7,15] 
[ 6,15] 


7.37/12 
15.96/14 


1.539+J 1 

1.3341?? 


[ 7,15] 
[ 7,15] 


23.62/12 
10.59/12 


6.0 


16 3 


•48 


0.13700 
0.13810 
0.13856 


0.817+^ 

0.678ig 7 

0.616+^ 


[10,23] 
[10,23] 
[10,23] 


21.26/22 
27.46/22 
26.95/22 


0.909±g 
0.810l^ 

0-774+21 


[10,23] 
[ 8,23] 
[ 8,23] 


23.49/22 
23.95/26 
41.53/26 


6.2 


24 3 


•48 


0.13640 
0.13710 
0.13745 


0.608+| 

0.509lg 2 

0.467l^ 


[11,23] 

[11,23] 
[11,23] 


34.87/20 
38.62/20 
22.68/20 


0.691+^ 

0.620tJo 
0.577±i§ 


[11,23] 
[11,23] 
[11,23] 


19.42/20 
23.87/20 
20.60/20 



TABLE VIII. Masses for S-like, A-like and A-like baryons for non-degenerate quark masses for 
(3 = 6.0, L 3 • T = 16 3 • 48 for the non-perturbatively improved data. 



«1 


K2 


^3 


arriY, 


am\ 


arriA 


0.13344 


0.13344 


0.13417 


0.780+^ u 


0.7751J 1 


0.894l^ 


0.13344 


0.13344 


0.13455 


0.766J3 1 


0.757l^ 2 


0.890±1| 


0.13344 


0.13417 


0.13417 


0.7431J 2 


0.744l^ 3 


0.8721J 4 


0.13344 
0.13344 


0.13417 
0.13455 


0.13455 
0.13455 


0.735±Jjj 

0.708+^ 


0.7221}^ 
0.71311? 


0.8711H 
0.860±1| 


0.13417 


0.13344 


0.13344 


0.771+J 1 


0.7811J 1 




0.13417 


0.13344 


0.13417 


0.7451* 4 


0.739lf 




0.13417 


0.13344 


0.13455 


0.7281} 2 , 


0.727lg 4 




0.13417 


0.13417 


0.13455 


0.6971JI 


0.692111 


0.845l|J 


0.13417 


0.13455 


0.13455 


0.677+]f 


0.679ll£ 


0.8371P 


0.13455 


0.13344 


0.13344 


0.750l^ 3 


0.769l| 2 




0.13455 


0.13344 


0.13417 


0.719±Jg 


0.7341} 3 




0.13455 


0.13344 


0.13455 


0.7171 J| 


0.7021} 4 




0.13455 


0.13417 


0.13417 


0.6891H 


0.698l^ 4 




0.13455 


0.13417 


0.13455 


0.6791JI 


0.674l^ 
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TABLE IX. Masses for S-like, A-like and A-like baryons for non-degenerate quark masses for 
j3 = 6.0, L 3 -T = 32 3 • 64 for the non-perturbatively improved data. 



«1 


K2 


^3 


arriY, 


arriA 


OlTlA 


0.13344 


0.13344 


0.13417 


0.769+S 


0.765l^ 


0.873l^ 


0.13344 


0.13344 


0.13455 


0.755+1° 


0.74811H 


0.8591H 


0.13344 


0.13417 


0.13417 


0.732+JH 


0.737+1°; 


0.8451H 


0.13344 


0.13417 


0.13455 


0.726+^ 


0.71111° 


0.8321H 


0.13344 


0.13455 


0.13455 


0.695+^ 


0.70411° 


0.820ll| 


0.13417 


0.13344 


0.13344 


0.764+?! 


0.77211? 




0.13417 


0.13344 


0.13417 


0.74011°, 


0.73011H 




0.13417 


0.13344 


0.13455 


0.7161 J2 


0.719l^ 




0.13417 


0.13417 


0.13455 


0.6841 J| 


0.680l^ 


0.8051H 


0.13417 


0.13455 


0.13455 


0.659lJ| 


0.6641H 


0.793l^ 


0.13455 


0.13344 


0.13344 


0.7451^ 


0.75911? 




0.13455 


0.13344 


0.13417 


0.712+1J 


0.724l^ 




0.13455 


0.13344 


0.13455 


0.7111}° 


0.6911JI 




0.13455 


0.13417 


0.13417 


0.678l^ 


0.68411? 




0.13455 


0.13417 


0.13455 


0.6661H 


o.657ll[; 





TABLE X. Masses for E-like, A-like and A-like baryons for non-degenerate quark masses for 
(3 = 6.2, L 3 -T = 24 3 • 48 for the non-perturbatively improved data. 



«1 


K2 


«3 


arriY, 


am\ 


arriA 


0.13460 


0.13460 


0.13510 


0.555+J 


0.5481? 


0.656+^ 


0.13460 


0.13460 


0.13530 


0.547ll 


0.534+J 


0.648+J 


0.13460 


0.13510 


0.13510 


0.528+^ 


0.524l| 


0.638l| 


0.13460 


0.13510 


0.13530 


0.52211°, 


0.510lJi 


0.630l^° 


0.13460 


0.13530 


0.13530 


0.513l[ 2 


0.502l? 2 


0.623+ 11 


0.13510 


0.13460 


0.13460 


0.545+J 


0.5521? 




0.13510 


0.13460 


0.13510 


0.525+J 


0.522l? 




0.13510 


0.13460 


0.13530 


0.514l? 


0.516l 10 




0.13510 


0.13510 


0.13530 


0.5001?! 


0.4971?! 


0.61lllg 


0.13510 


0.13530 


0.13530 


0.49711! 


0.474l^ 7 


0.60611? 


0.13530 


0.13460 


0.13460 


0.5301? 


0.543l^° 




0.13530 


0.13460 


0.13510 


0.50911! 


0.51611J 




0.13530 


0.13460 


0.13530 


0.5031?! 


0.507l? 4 




0.13530 


0.13510 


0.13510 


0.493ll 2 


0.497113 




0.13530 


0.13510 


0.13530 


0.487ll 3 


0.497+f 5 
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TABLE XL Masses for S-like, A-like and A-like baryons for non-degenerate quark masses for 
(3 = 6.0, L 3 -T = 16 3 • 48 for the tadpole improved data. 



«1 


K2 


^3 


amj 


am^ 


OlTlA 


0.13700 


0.13700 


0.13810 


0.777+^ 


0.7711J 1 


0.873l^ 


0.13700 


0.13700 


0.13856 


0.761^ 2 


0.7491J 3 


0.882lf 5 


0.13700 


0.13810 


0.13810 


0.72213, 4 


0.732l^ 3 


0.853l 4 8 


0.13700 


0.13810 


0.13856 


0.713^ 5 


0.703l^ 6 


0.845l? 7 


0.13700 
0.13810 


0.13856 
0.13700 


0.13856 
0.13700 


0.678iijj 

0.76615 1 


0.6971J 5 

O.778I5 2 


0.832l{^ 


0.13810 


0.13700 


0.13810 


0.733l£ 4 


0.718l^ 4 




0.13810 


0.13700 


0.13856 


0.706l^ 4 


0.707l^ 5 




0.13810 


0.13810 


0.13856 


0.6591 J| 


0.659l 2 ^ 


0.8101^ 


0.13810 
0.13856 


0.13856 
0.13700 


0.13856 
0.13700 


0.6341^2 
0.7411J, 3 


0.643l 2 ^ 
0.7631J 4 


0.787111 


0.13856 


0.13700 


0.13810 


0.698l^ 7 


0.710l^ 6 




0.13856 


0.13700 


0.13856 


0.697l^ 5 


0.6721H 




0.13856 


0.13810 


0.13810 


0.6541 2 ] 


0.6551J1 




0.13856 


0.13810 


0.13856 


0M2±H 


0.626l^ 





TABLE XII. Masses for S-like, A-like and A-like baryons for non-degenerate quark masses for 
(3 = 6.2, L 3 • T = 24 3 • 48 for the tadpole improved data. 



«1 


K2 


«3 


arriY, 


arriA 


arriA 


0.13640 


0.13640 


0.13710 


0.576!^ 


0.5681^ 


0.656l£ 


0.13640 


0.13640 


0.13745 


0.5541° 


0.545l^ 


0.6521? 


0.13640 


0.13710 


0.13710 


0.5391^ 


0.5351J 1 


0.6421? 


0.13640 


0.13710 


0.13745 


0.526l^ 


0.510ljj 


0.6271 1 , 


0.13640 


0.13745 


0.13745 


0.504l5 , 3 


0.491±lg 


0.6191J 1 


0.13710 


0.13640 


0.13640 


0.5631!? 


0.5711? 




0.13710 


0.13640 


0.13710 


0.539 + £ 


0.5321J 1 




0.13710 


0.13640 


0.13745 


0.515l^° 


0.518l| 




0.13710 


0.13710 


0.13745 


0.494l^ 


0.4911^ 


0.60211J 


0.13710 


0.13745 


0.13745 


0.463l^ 


0.457l| 2 


0.59311? 


0.13745 


0.13640 


0.13640 


0.540+^ 


0.5581? 




0.13745 


0.13640 


0.13710 


0.5101*° 


0.5221*° 




0.13745 


0.13640 


0.13745 


O.497I9 5 


0.492l? 5 




0.13745 


0.13710 


0.13710 


0.4701^ 


0.4921^3 




0.13745 


0.13710 


0.13745 


0.4551^3 


0.45ll| 
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TABLE XIII. Values for n c for all datasets. The labels 'tree', 'bare' and 'tad' refer to the 
coefficient b m estimated either at tree-level, bare and tadpole improved one-loop perturbation 
theory, respectively. Also shown are the results of ref. ||, obtained using the current quark mass 
instead of the pseudoscalar mass. 








L 3 


T 


6m = 


tree 


bare 


tad 


ref. § 


NP 


6.0 
6.0 


16 3 
32 3 


48 
64 




0.135259lf 

0.13524ll? 


0.135255i£ e 

0.135237l? 


0.1352521^ 
0.135235l? 


0.135196(14) 




6.2 


24 3 


48 




0.135818+^ 


0.135816+^ 


0.135815+^ 


0.135795(13) 


TAD 


5.7 


16 3 


48 


0.143408+^ 


0.143240±jjg 


0.143206+^ 


0.143179ig| 






6.0 


16 3 


48 


0.139240+7° 


0.139216^ 9 


0.139212+t 9 


0.1392091^ 






6.2 


24 3 


48 


0.137912±];| 


0.137900+^ 


0. 137898+ J| 


0.1378971^ 





TABLE XIV. Values for n n and k s determined for three different quantities Q to set the lattice 
scale. All estimates were obtained by setting b m equal to its tadpole improved perturbative value. 



Q 



Q = rriK* 



Q = m N 



13 L 3 



K 7i 



0.13520+ 1 0.13401+^ 



0.133181^ 
0.13327+.$ 
0.134381 1 ! 



NP 



6.0 
6.0 
6.2 



16 3 
32 3 

24 3 



IS 



0.13520] 



64 0.13519+] 0.13398+J 0.13518 
48 0.13578+ 1 0.13495+J 



48 0.14306+^ 0.1401211 



i 

i 

-i 

0.135771? 



0.14307lf 

0.139151? 

0.137851? 



0.133831^ 
0.133761^ 
0.134761 3 , 



0.140291^" 

0.137691^ 

0.13680 



0.1351715 

0.13516l{ 
0.135761? 



0.14302+ 3 
0.139131? 
0.137841? 



0.139041^? 
0.13712+ 13 , 
0.13655+ 1 ? 



TAD 



5.7 
6.0 
6.2 



16 3 
16 3 

24 3 



48 



0.139151^ 0.13780+ 1 



48 0.13786 



-2 



0.13699 



+2 
-2 



-5 
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TABLE XV. Results for the physical meson and baryon masses with the lattice scale set by ro- 



m p r Q 
(3 NP TAD 



m K *r 
NP TAD 



NP TAD 



-54 



2.224+J* 
2.422j^ 2.356lf| 



2.4661^ 
2.655±|| 2.592±|| 

2.682±|| 2.654JH 

2.729(77) 



5.7 
6.0 

6.2 
Cont. 



1.982; 



2.189±|| 2.12lj^ 

2.254+1? 2.211±$ 

2.352(163) 



2.468 



+51 
-55 



2.433 



+48 
-50 



2.540(117) 



m N r 
(3 NP TAD 



msro 
NP TAD 



m A r 
NP TAD 



771= Tq 

NP TAD 



5.7 

6.0 

6.2 

Cont. 



3.08] 
3.02^ 



-13 
-13 



11 



7 



2.74 H 



2.92+^ 



2.87 H 



12 



3.36 
3.31" 



-10 

-li 



3.05 
3.22 



ii 



- 5 

11 

- 7 

3.181 I 





3.05+ I 


3M±\\ 


3.21+ 1 * 


3.30^ 


3-lSl/o 



2.92(24) 



3.23(19) 



3.22(20) 



3.63] 



9 
8 

3.59+ I 



3.351 
3.52" 



3.501 7 



3.54(15) 



mAro 
NP TAD 



m E *r 
NP TAD 



777=* Tg 

NP TAD 



3.29l^ 

-ii 



3^5+T 
4.131J 

9 

7 



""iO^O 

NP TAD 



5.7 

6.0 

6.2 

Cont. 



4.12 



-26 
-19 



3.93: 



■02+g 



16 



3.96!^ 



-20 
-16 



3.86(37) 



4.29; 

4.241$ 4.18 
4.15(29) 



3.811 g 

4.47l^ 4.32l^ 

4.46l 1 g 4.401 e 

4.44(22) 



4.641 I 



4.081 I 



4.52: 
4.681 6 4.62^ 
4.72(17) 



TABLE XVI. Same as Table XV but with the scale set by tuk 



d 



m p /m K * 
NP TAD 



NP TAD 



5.7 
6.0 
6.2 

Cont. 



a 



0.8961 9 
0.8901$ 0.8871$ 

903+ 12 n snQ "'" ] 



-25 



0-92111 



0.893133 



1.1131 7 
1-H2±i§ 



1.105; 

1.107; 
1.107; 



- 2 

- 4 

- 4 

- 6 

- 2 
-10 



1.110: 
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m N /m K * 
NP TAD 



NP TAD 



m A /m K * 
NP TAD 



NP TAD 



^3T~ 
-37 
+57 

-41 



+37" 

-29 

-16 

-32 

1-3121JI 



1.372l^ 
1.356±£? 

1.313l^ 



T504+1 
1.5001H 

1.465l| 



5.7 

6.0 

6.2 

Cont. 



/3 



1.239 
1.220 



1.253+Ji 

1.212±i| 1.16011 



l.«385_gn 

1.35211 



1.372 
1.360 



1.14+ 



18 



1.291^ 



1.347li 



1.29+J 



1.5171^ 

1.4911^ 



1.45!^ 



NP TAD 



NP TAD 



771s* /m^* 

NP TAD 



m n /m K * 
NP TAD 



1.4871^ 
1.672199 1.645l^ 



1.600l^ 
1.7561$ 1-7351^ 



I^Islf" 

1.84ll^ 1.825+J5 
1.823±li I.8HI32 

1 79+ 9 



^3TT 
•21 



5.7 

6.0 

6.2 

Cont. 



1.826' 



1.60911 



I.598I47 



1.501J? 



l-716li 
1.64 



1.705l| 

+13 



1.926l|§ 1.916lj$ 



-13 



1.93011 
1.93 



1.918118 

+ 7 
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TABLE XVII. Same as Table |XV| but with the scale set by the nucleoli mass tojv- 





m p /m N 


rriK* /tun 


m^/m N 




P 


NP TAD 


NP TAD 


NP TAD 




5.7 


0.727±3Q 


0.8431^ 


0.959l}^ 




6.0 


0.709JH 0.726±§J 


0.837l^ 0.848l|? 


0.9651$ 0.970l^ 




6.2 


0.744+^ 0.769+J? 


0.859+ft 0.882l^ 


0.975l 2 7 1 0.995l^ 




Cont. 


0.796+j| 


0.894+Jg' 


0.992l|^ 






ms/m N 


m\/m N 


m= /mjv 




P 


NP TAD 


NP TAD 


NP TAD 




5.7 


1.1501 I 


1.1501 3 


1.2991 5 




6.0 


1.151+lg 1.1561 ? 


1.148± 1 § 1.1531 g 


1.303l$j 1.31lli3 




6.2 


1.1541 12 1.1621 | 


1.15ll^ 1.1621$ 


1.308lf| 1.3231 1 ! 




Cont. 


1.16111 


1.157115 


1.321+J1 






m^/uiN 


mY,*/m N 


7n=* /mjv 


m n /m N 


P 


NP TAD 


NP TAD 


NP TAD 


NP TAD 


5.7 


1.198tg 


l-326±$ 


1.4541* 


1.582l^ 


6.0 


1.330±?§ 1.344lg 


1.4271$ 1.445l^ 


1.524l|? 1.545l§| 


1.62ll$ 1.6461H 


6.2 


1.3251^ 1.37611 


1.44311 1.488131 


1.561+lg 1.601+i 


1.679+^ 1.714 + |° 


Cont. 


1.33±?g 


1.471 1 ? 


1.611 1 ! 


1.761^ 



TABLE XVIII. Results for the physical meson and baryons masses at (3 = 6.0 on 32 3 ■ 64. 
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TABLE XIX. Results for vector mesons, octet (S-like) and decuplet (A-like) baryons, in units 
of 7"o, interpolated to the reference points defined by (mpsro) 2 = 3.0 and mps/rn\r = 0.7. 
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FIG. 1. The vector-pseudoscalar hyperfine splittings for (a): theNPand(b): the TAD datasets. 
Open squares, full circles and crosses denote the data at j3 = 6.2, 6.0 and 5.7, respectively. The 
experimental points are represented by the asterisks. 
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FIG. 2. The parameter J plotted versus the lattice spacing in units of ro- Filled (open) symbols 
denote the data using non-perturbative (tadpole improved) c sw . The filled square denotes the data 
point at j3 = 6.0 on the larger volume of 32 3 • 64. 
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FIG. 3. The pseudo-scalar squared plotted versus the averaged quark mass, m q = ifht^x+fh^z) /2 
in units of vq for (a): the NP and (b): the TAD datasets. Open squares, full circles and crosses 
denote the data at (3 = 6.2, 6.0 and 5.7, respectively. 
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FIG. 4. Data for (a): vector mesons, (b): X-like baryons and (c): A-like baryons plotted versus 
the averaged quark mass m„ = (ifia,i + fhq,2 + ^q,3)/3 in units of ro for the non-perturbatively im- 
proved (NP) dataset. Open squares and full circles denote the data at (3 = 6.2 and 6.0, respectively. 
The lines represent the fits to eqs. (p^-18). 
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FIG. 5. Same as Fig. || for the TAD dataset. Open squares, full circles and crosses denote the 
data at j3 = 6.2, 6.0 and 5.7, respectively. 
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FIG. 6. Examples of simultaneous continuum extrapolations of the NP and TAD datasets. (a): 
the K* meson, (b): the nucleon and (c): the A. Full circles denote the data computed using 
the non-perturbative estimate of c sw , whereas the tadpole improved data are represented by open 
squares. 
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FIG. 7. The quenched light hadron spectrum computed in the 0(a) improved and its compari- 
son to the results of ref. Q, obtained using the unimproved Wilson action (full circles). The levels 
of the experimental points are denoted by the solid lines. 
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